Partially Asymmetric Exclusion Processes with Sitewise Disorder 
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We study the stationary properties as well as the non-stationary dynamics of the one-dimensional 
partially asymmetric exclusion process with position dependent random hop rates. In a finite system 
of L sites the stationary current, J, is determined by the largest barrier and the corresponding 
waiting time, r ~ J , is related to the waiting time of a single random walker, 
The current is found to vanish as: J ~ L~ Z//2 , where z is the dynamical exponent of the biased single 
particle Sinai walk. Typical stationary states are phase separated: At the largest barrier almost all 
particles queue at one side and almost all holes are at the other side. The high-density (low-density) 
region, is divided into ~ L 1 ^ 2 connected parts of particles (holes) which are separated by islands of 
holes (particles) located at the subleading barriers (valleys). We also study non-stationary processes 
of the system, like coarsening and invasion. Finally we discuss some related models, where particles 
of larger size or multiple occupation of lattice sites is considered. 



I. INTRODUCTION 

The stochastic dynamics of self-driven particles on one- 
dimensional lattices is one of the key problems of non- 
equilibrium physics 1 ^. Models of this kind have been 
used in order to describe such different problems as e.g. 
highway traffic^ or the dynamics of motor proteins along 
actin filaments or microtubuli 4 . These two examples have 
in common that the directed particle motion against dis- 
sipative forces is maintained by the steady input of en- 
ergy, a feature that leads to the generic non-equilibrium 
behaviour of the system. The finite drift of the particles 
leads to a strong sensitivity to spatial inhomogeneities. 
In contrast to equilibrium systems, even a single local 
defect may lead to bulk effects as e.g. the separation 
into macroscopic low and high density domains^. More- 
over disorder of any type may also strongly influence the 
transport capacities of the system^. This is in particu- 
lar true for the case of strong disorder, i.e. realizations of 
the disorder where the local direction of the bias is non- 
uniforn&^iiSiii. The prototype of this kind of stochas- 
tic motion is the Sinai walk, which has been studied to 
great extencUi. The stochastic motion of a Sinai walker 
is characterized by large velocity fluctuations. The ori- 
gin of these fluctuations is most easily understood if one 
translates the local hop rates into energy differences, i.e. 
a forward bias corresponds to a negative slope of the cor- 
responding energy landscape and vice versa. Obviously 
the energy landscape is not flat for strong disorder, but 
characterized by barriers of different heights. The time a 
walker spends in front of such a barrier increases expo- 
nentially with the height of the barrier and therefore the 
largest barriers determine the behavior of the system. 

There exists a large class of both quantum and stochas- 
tic models where the disorder plays a dominant role over 
deterministic fluctuations and where a connection to the 
Sinai walk can be made 1 ^. In this respect one can men- 



tion random quantum spin chains 1 ^* 1 ^ or reaction diffu- 
sion models with quenched disorder—. Besides a number 
of experimental setups have been identified, where this 
model properly describes the particle dynamics, e.g. the 
translocation of a RNA strand through a porei& or the 
motion of molecular motors on microtubules^. 

While most of the models refer to the single parti- 
cle case, more recently, the effect of strong disorder on 
driven-diffusive many particle systems, i.e. the so-called 
asymmetric exclusion process with strong disorder, has 
been investigatec&SiiSiiLiSiiSiSii. It has been shown that 
the properties of the system largely depend on the way 
how the disorder is implemented. In case of particle dis- 
order, i.e. for the case where the hopping rates of the 
particles are quenched random variables—, the transport 
properties of the system are in close analogy to the sin- 
gle particle systerr&iSi both in the Griffiths phase, i.e. 
in case of a tilted energy landscape, and at the critical 
point, i.e. for a vanishing drift velocity of the particles. 
This implies that the many particle effects are most pro- 
nounced for quantities that are related to the arrange- 
ment of the particles. In case of lattice disorder, however, 
the transport capacity of single and many particle sys- 
tems are also differenljSiS, because the presence of many 
particles alters the effective heights of the energy barri- 
ers. A relation between the stationary current for particle 
disorder and for lattice disorder is announced in Ref.0. 
In the present work we elaborate the properties of the 
distribution of the stationary current and show connec- 
tions with extreme value statistics2i*2£. We also discuss 
the properties of the density profile as well as the conse- 
quences for the coarsening behavior. Some of the results 
obtained in this work may also be relevant for models 
which are closely related to the exclusion process, such 
as the Heisenberg chainSi and the dimer evaporation and 
deposition process^. 

The article is organized as follows. The model is in- 
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troduced in Sec ITU and the basic properties of the single 
particle motion in the framework of the concept of trap- 
ping or waiting times is given in Sec lIIII The stationary 
properties of the many particle motion is prestented for 
non-zero average bias in Sec II VI and for zero average bias 
in SecE] Some related models are discussed in Sec I VII 
and non-stationary phenomena is described in Sec lVIII 
In the final section of this article we will summarize and 
discuss our results.. 



II. MODEL 

We consider the partially asymmetric exclusion process 
(PASEP) with site-dependend hop rates on a one dimen- 
sional lattice. Each lattice site i can either be empty, i.e. 
n = 0, or occupied by a single particle, r» = 1. 

The time evolution of the local configuration (Ti,Tj+i) 
is described by 
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with rate Pi, 
with rate 



(1) 

(2) 



i.e. particles are hopping in positive direction (from i to 
site i + 1) with rate pi and in negative direction (from 
site i + 1 to site i) with rate qt. 

It is useful to relate the stochastic dynamics to a po- 
tential landscape by the relation 



* = e -{Ui-U i+ r) 



(3) 



where t/, is the energy assigned to site i (relative to a 
reference value, U\ = 0). A typical sample of this type 
of landscape is shown in Fig. [21 Links with backward 
bias (qt > Pi) correspond to a descending links of the en- 
ergy landscape while those with forward bias to ascend- 
ing ones. It turns out that the large scale behavior of 
the process is determined by the wandering properties of 
the energy landscape. The relation between energy land- 
scape and system properties will be discussed in some 
detail in the next section. 

The hop rates are independent and identically dis- 
tributed random variables taken from the distributions, 
p(p)dp and Tr(q)dq, respectively which will be specified 
later. We restrict ourselves to types of randomness where 
there are both forward and backward links with finite 
probability, i.e. the easy direction of hopping is also ran- 
dom. We define a control parameter as the average asym- 
metry between forward and backward rates as follows: 



[Inp] 

av 

[Inq] 

av 

var[lnp] + var[ln q] 



(4) 



such that for S > (5 < 0) the particles move on av- 
erage to the right (left). Here, and in the following 
[. . . ] av denotes average over quenched disorder, whereas 
var(z) stands for the variance of x. Note, that S ~ 
limi—too Ul/L, thus the control-parameter is propor- 
tional to the average slope of the potential landscape. 



In most of the calculations we use a binary disorder 
distribution: 

p(p) = c6(p — 1) + (1 — c)5(p — r) 

Piqi = T < 1 for all i , (5) 

where the control-parameter takes the value: 

1 - 2c 



5 = 



2c(l — c) lnr ' 



(6) 



thus a biased motion to the right is realized for 1 < c < 
1/2. 

In some cases we use a power-law distribution: 

P(P) = -jjP~ 1+1/D . for0<p<l; (7) 



-l/D 



where D 2 — var[lnp] = var[lng] measures the strength of 
disorder and S = \nq / (2D 2 ). 



III. SINGLE PARTICLE MOTION - THE 
CONCEPT OF TRAPPING TIMES 

We wish to introduce in the next section a simple phe- 
nomenological approach to the the biased exclusion pro- 
cess, which is based on the dynamics of the correspond- 
ing single-particle problem, known as Sinai walkii^. It 
is a thoroughly studied problem and here we recapitulate 
some of its properties which are necessary to understand 
to many-particle problem. For a periodic chain of length, 
L, the stationary drift velocity, vq, for a given realization 
of disorder is given by2£: 



1 

Vq 



(8) 



where 



^=i i+ En— 



(9) 



is the persistence probability^ at site This quantity 
measures the fraction of walks which start at site i + 1 
and pass the i, i + 1 link first from the left, i.e. after a 
complete tour along the chainSi. 

It is clear that the regions where the particle spends 
long times are valleys in the energy landscape which are 
followed by a large barrier represented by a Brownian 
excursion, see the section is) in Fig. ^ It is obvious 
that one can define a (single particle) trapping time for 
each barrier (Brownian excursion): 
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(10) 
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FIG. 1: Schematic energy landscape. 



which is the time needed to escape from that barrier and 
where i is the starting point of the corresponding excur- 
sion, i.e. the minimal energy site of the valley in front 
of the barrier (site i\ in Fig. QJ. The trapping time is 
of the form of a Kesten random variable 2 ^, moreover one 
can convince oneself that trapping times belonging to dif- 
ferent barriers are practically independent. It is known 
that in the large- L limit the probability distribution of 
such variables has an algebraic tail: 



P(t) 



.-1-1/2 



(11) 



where the exponent, z, is the positive root of the equation 



l/z 



= 1. 



(12) 



Choosing the bimodal disorder as defined in Eq.JSJ) one 
obtains: 



lnr 



ln(c-! - 1) 



(13) 



Now concerning the drift velocity of the random walker 
it is the inverse of the average trapping time for z < 1. If, 
however, z > 1, the average of the trapping time is diver- 
gent and the motion of a single particle is determined by 
the largest trapping time, r m , the typical value of which 
follows from the relation^: / T °° P(t)AtL = 0(1) and 
given by r m ~ L z . 

The persistence probability and thus the trapping time 
at a given barrier is related to the potential landscape and 
for large trapping times it is asymptotically given by: 



L-l 

E 

fe=i 



i 



Pi+k 



-Ui 



(14) 



which is well approximated by the largest term in the 
sum: Ui = max fc ([/i +fc - Ui). 

In this way the motion of the particle for z > 1 is 
influenced only by the deep valleys and the correspond- 
ing large barriers, and ultimately its velocity is deter- 
mined by passing the largest barrier in the system, since 



~ e 



and U„ 



max?/,. The intimate relation be- 



of the energy landscape led to the idea of renormaliz- 
ing the energy landscape 12 '—. By using the renormaliza- 
tion scheme and alternative approaches one obtains the 
same form of the trapping time distribution 30 as given in 
Eq. (TTJ. 



IV. MANY PARTICLE MOTION 

We now turn to the many particle case, i.e. we put 
N = 0(L) particles on the lattice. First we consider a 
special form of disorder for which exact results can be de- 
rived. These are then generalized for more general form 
of disorder using phenomenological and scaling consider- 
ations. 



A. Extreme binary disorder 

Here we consider the bimodal disorder as defined in 
Eq.© where the landscape is represented by a random 
walk having a step of size lnr with probability (1 — c) 
and of size — In r with probability (c) . In the limit c <C 1 
the walk is strongly biased, having typically downwards 
steps and the rare upward steps are forming the barri- 
ers. If we also have r -C 1 the large barriers are typically 
straight, since for a large energy-scale the fluctuations 
are reduced—. The largest barrier consists of l m upward 
steps, has a height of U m = l m In r _1 and we have for the 
typical value c lm L — 1. It is easy to check that the trap- 
ping time for single particle motion is r m ~ r lm ~ L z 
and the dynamical exponent z' = lnr/ In c -1 corresponds 
to that in Ea. ll3|) in the given limit. Also for many par- 




tween the mobility of the particle and structure of the 



FIG. 2: Schematic form of the potential landscape for ex- 
treme disorder and the typical position of the particles. The 
high-density and low-density phases are separated at the mid- 
dle of the largest barrier (2), the arrow indicates the height 
Um/2. Deep valleys in the low-density phase are filled un- 
til the effective barrier is Um/2 (3) and high barriers in the 
high-density phase are vacant above a level of C/ m /2 (1), see 
text. 

tides the largest barrier plays the dominant role and we 
analyze first the stationary current in a system having 
just one large barrier of length l m . This problem corre- 
sponds to a PASEP with particle in- and output against 
the direction of the bias. In the stationary state of this 
model, which is known exactly^!, the system is half-filled 
and there is a sharp front in the middle of the chain, 
which is the consequence of particle-hole symmetry, so 
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that particles (and holes) should overcome only an ef- 
fective (filled) barrier of height U m /2. Consequently the 

1 /2 

effective trapping time is f m ~ r m and the stationary 
current is given by: 

J~L- Z/2 (15) 

which goes to zero in the thermodynamic limit. 

Now, if we consider the complete potential landscape 
the basic features of the above considerations remain 
true, since the largest barrier governs the stationary dy- 
namics of the system. In particular in the stationary 
state there is a sharp front located at the middle of 
the largest barrier. Thus, due to particle-hole symme- 
try the largest barrier is half-filled and the stationary 
current is related to the largest effective trapping time 
and given by Ea. (|15|l . The subleading barriers lead to 
a substructure both in the high-density and in the low- 
density phases, which can be related to each other by 
using the particle-hole symmetry. In the low-density re- 
gion it is easy to see that small barriers with Ui < U m /2 
are completely empty, since they do not cause sufficient 
resistivity for the current. Only sufficiently large barriers 
with Ui > U m /2 are able to slow down the particle and 
cause jams. These barriers will be filled up to the level: 
Ui — U m /2 as a consequence of the conservation of the 
stationary current. The typical number of occupied bar- 
riers in the low-density phase is given by: c lm ^ 2 L ~ L 1 ' 2 . 
In the high-density phase due to particle-hole symmetry 
holes are accumulated only in such subleading barriers 
for which Ui > U m /2. The number of holes in such a 
barrier scales as Ui — U m /2 and the typical number of 
occupied barriers is O^ 1 / 2 ). 



B. General form of disorder 

Next we extend our discussion to general form of dis- 
order, when it may be still assumed that there is a sharp 
front separating the high-density and low-density regions. 
If this is the case, the front is located at a position, where 
the escape rate of particles and the escape rates of holes 
which leave the front to the left through the high-density 
region, are equal, which is a necessary condition for the 
the localisation of the front. As in the extreme disorder 
case, the front is located in the middle of the largest bar- 
rier.. The potential landscape compared with that for 
extreme disorder in Figgis modified in such a way that 
instead of long straight lines the barriers and valleys are 
made by Brownian excursions. 



1. Distribution of the current 

The typical value of the stationary current follows from 
the same reasoning as in Section llVAl it is given by the 
square root of the trapping time of a single walker at 



the largest barrier. Furthermore the trapping time dis- 
tribution follows an asymptotically algebraic distribution 
given by ff 1 1 11 , 

Therefore for a given system we have to find the maxi- 
mum out of 0{L) algebraically distributed random vari- 
ables, a problem which is thoroughly studied in the math- 
ematical literature (see e.g. in |21|) and has demonstrated 
to work for strong Griffiths singularities in disordered 
systems 2 ^. For our case it follows that the current is de- 
scribed by the well-known Frechet distribution given by 

P(j) = lp/*-i e -3 2/ \ (i6) 
z 

in terms of J = JqJL z / 2 , where the non-universal con- 
stant Jo depends on the prefactor of the tail. Thus, the 
current scales with the system size as given in Eq. (|15|l . 

Although the above results hold for generic hop rate 
distributions, a caveat is in order concerning the type of 
randomness. Besides the bimodal one, one may also con- 
sider hop rate distributions where the support of the the 
forward rates, p, does not have a positive lower bound, 
such as for the power-law distribution in Eq.Q. In this 
case it may happen that the current is limited by the 
lowest hop rate p m i n ~ L~ D , i.e. by a single slow link 
rather than by the maximal amplitude of a Brownian ex- 
cursion, which is given in Ea. H15l) . Clearly, the exponent 
which enters in (|16|l and as well as in the scaling relation 
of the current (|15f) is now max{z/2, D}. 

For a uniform distribution with D = 1 the current 
scales as J(L) ~ L~ 2 / 2 for z > 2 whereas J(L) ~ L _1 
for 1 < z < 2. Numerical results for the distribution of 
current in the anomalous region 1 < z < 2 are shown in 
Fig. whereas numerical results for the bimodal ran- 
domness, for which the anomalous scenario never sets in, 
can be found in Refit In what follows we assume that 
z/2 > D always holds and the current is controlled by an 
extended region. 




-2 2 4 6 

x=-ln(JL) 

FIG. 3: The distribution of the current, calculated with 
uniform randomness with qo = 1/3, where z — 1.335. The 
number of samples is 10000. The solid line is the Frechet 
distribution given in 1161 . 
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2. Density profile 

The density profile for general form of disorder has the 
same qualitative features as for the extreme disorder in 
Sec. IIV Al The system consists of a high-density phase, 
where all the lattice sites are typically occupied, and of 
a low-density phase in which all the lattice sites are typi- 
cally vacant. The front separating the two phases is just 
at the middle of the largest barrier. In the low-density 
phase only the large valleys with Ui > U m /2 are filled to 
a level Ui — U m /2. Similarly, in the high-density phase 
only the large barriers with Ui > U m /2 are only partially 
filled up to U m J 2. 

The number of occupied valleys n(L), (and the number 
of partially filled barriers) can be estimated from the con- 
dition, that at these barriers the single particle trapping 
time is larger than 1/J, thus 

/>oo 

n{L) ~L p{r)dT ~ LJ 1 / 2 , (17) 
Ji/J 

where we used Ijlljl and the fact that, due to their finite 
typical extension, the number of barriers in the system 
is of order L. Keeping in mind the scaling relation of 
the current in (|15|l . we obtain n(L) ~ L 1 / 2 , in agreement 
with the extreme disorder case in Sec. IIV Al This finding 
is in good agreement with numerical results shown in Fig. 

El 

The length of Brownian excursions has an exponen- 
tial distribution 32 , with a finite ^-dependent characteris- 
tic value. Although, the length of the L 1 / 2 longest ex- 
cursions is 0(lnL)2i, the length of the empty domain in 
each of these barriers equal approximately to the half of 
the length of the largest barrier, so one can show that 
the number of particles they contain is typically finite. 
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FIG. 4: Numerically calculated stationary density profile of 
a sample of size L — 512. The hop rates were generated with 
the bimodal distribution with c = 0.3 and r = 0.5 in Eq.@. 

The total number of particles clustering behind the 
subleading barriers is thus of the order of L 1 ^ 2 . The 
other O(L) particles must be behind the largest barrier, 
where a macroscopic cluster of particles is formed. Ac- 
cording to the discussion given above, the particles can 




FIG. 5: Average number of points where the density profile 
crosses the line p(x) — 1/2 calculated numerically for different 
system sizes L. Binary randomness were used with r = 0.5 
and c = 0.2, whereas the average was performed over 100 
samples. The straight line has slope of 1/2. 

be categorized into three classes as was carried out in a 
similar way for the disordered zero-range processes. 

Almost all particles are found at the largest barrier, 
forming a condensate, where the density is close to one. 
The length of the condensate and the number of particles 
in it is thus 

A^ c ~ N. (18) 

This region is, however, cut into sections of size ~ L 1 / 2 
by vacant clusters, which have a finite length. 

In the remaining part of the system which is of size 
~ [L — N), particle clusters of finite size are found, sep- 
arated by a typical distance of L 1 / 2 . Due to the small 
local velocity in these dense clusters, the contribution 
of a single particle to the current is small compared to 
freely moving particles, therefore we call these particles 
inactive. Their total number is given by 

N ia ~ L 1 ' 2 . (19) 

The third class of particles, i.e. active ones, move 
practically freely between clusters of inactive particles. 
Consequently, they have a significant contribution to the 
current. In order to estimate the number of active par- 
ticles, we notice that due to the low particle density 
within sections between two neighboring clusters of in- 
active particles the dynamics of active particles corre- 
sponds to the dynamics of a single random walker in a 
disordered landscape. The typical length of these sec- 
tion is ~ L 1 / 2 and the largest trapping time of an active 
particle is ~ ^pj\ ~ L z / 2 . Now if z < 1, the particle 
velocity v is L-independent. Making use of the finite- 
size scaling form of the current, we obtain for the density 
p = J/v ~ L~ z / 2 , which is indeed vanishing in the large 
L limit. The total number of active particles is thus 

N a ~ L 1 -^ 2 , (z < 1). (20) 

For z > 1, the largest barrier in the empty section de- 
termines the travel time and the average velocity of par- 
ticles scales according to v(L) ~ L^~ z ^ 2 . The density 
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thus scales as p = J/v ~ L~ x / 2 , which means that typi- 
cally 0(1) particles reside between two inactive clusters. 
In this case the total number of active particles is given 
by 

N a ~ L 1 ' 2 , (z > 1). (21) 

According to the discussion above, the distribution of 
interparticle distances in a finite system of size L has a 
finite cut-off: they have the scaling form J ~ L z / 2 for 
z < 1, whereas I ~ L 1 / 2 for z > 1. 

Next, the number of particles at the second largest bar- 
rier ri2 is analyzed. The height of the 2nd largest barrier 
is proportional to InL just as the height of the largest 
one. For the bimodal randomness where the increments 
of the potential are fixed, the length of the barrier must 
grow with the system size at least as OQnL). On the 
other hand, it is generally true for types of randomness 
under consideration, where the current is controlled by 
an extended cluster of exponentially rare links. The typ- 
ical value for (712) is thus OQnL). In certain samples, 
however, where the single-particle trapping times T\ and 
t 2 at the 1st and 2nd largest barriers, respectively, are 
almost degenerate, meaning that r 2 < t\, the occupa- 
tion (n 2 ) is much larger, than the typical value O(lnL). 
We show below, that these samples give a singular con- 
tribution to the sample-averaged occupation [(n,2}] a v In 
samples, where the two largest barriers are almost degen- 
erate, n<i has large-scale fluctuations and the dominant 
contribution to the thermal average (n.2) comes from the 
fluctuations where n-i is much larger than its most prob- 
able value. In the latter case half-filling is realized at the 
2nd barrier, leading to an escape time ^/tJ. Not consider- 
ing the improbable samples with higher degeneracies (the 
contribution of which has in fact the same type of singu- 
larity), the number of particles located at barriers other 
than the two largest ones is negligible, and the problem is 
reduces to an effective two-barrier process with hop rates 
yj~¥\ and s/ri- Taking into account, that the total num- 
ber of particles is N, we obtain for the thermal average^ 



(n 2 )(a) = - (N + l)- TyTT; a<l, (22) 



where a = y/ ti jr\ , and (na) = N/2 for a = 1. Using 
the distribution function, p{a) we can average over real- 
izations: [(«2)]av = Jq { n )( a )p{ a )da, which is dominated 
by the contribution as a — > 1, where p{a) has a finite 
limiting value. Keeping in mind that the maximal value 
of (77,2) is N/2, we can write 

<n 2 ) av « / p(a)da~ p(l)hxN. (23) 

JO \-OL 

Thus the rare samples, where the second largest trapping 
time can be arbitrarily close to the largest one, lead to 
a logarithmically diverging contribution to the sample- 
averaged occupation at the 2nd barrier. 



C. Low density limit 

Contrary to the previous section, where the global den- 
sity p = N/L was finite, we study here the case, where 
the number of particles scales as N ~ CL a for large L, 
with < a < 1. This implies that the particle density 
is vanishing algebraically according to p(L) ~ CL a ~ x as 
L — > 00. In this case the number of active particles may 
be limited by N, which influences the current. 

We have learned in the previous chapter, that the sys- 
tem has a limited capacity for storing particles, given 
by N a + Ni a , and the excess of particles is driven to the 
macroscopic condensate. Evidently, as long as N exceeds 
the capacity of the system, N a + Ni a , the excess particles 
accumulate in the condensate and the largest barrier is 
haf filled. Consequently the current remains the same as 
the one observed for constant density. This is, however, 
no longer true if N is smaller then the capacity of the 
system. 

For z > 1, the capacity scales as (N a + N ia ) ~ L 1 / 2 
according to J3lJ and ifTi^l . Thus, for a > 1/2 the cur- 
rent, the number of active and inactive particles are not 
influenced, whereas N c ~ L a . However, if a < 1/2, the 
condensate vanishes, and since one can show that N a is 
proportional to N ia for arbitrary current, both must be 
proportional to L a . Therefore according to (|17|) the cur- 
rent is given by J - (n(L)/L) z - L _2(1 - a) . 

For z < 1, the capacity is dominated by N a , which 
scales according to (|2U|) . Thus, as far as a > 1 — z/2, the 
number of active and inactive particles and the current 
agrees with the case of a finite density, whereas the size 
of the condensate scales as JV C ~ L". If a < 1 — z/2, 
there is no condensate present in the system and the 
number of active particles is limited to N a ~ L a . There- 
fore the corresponding current scales as J ~ L~( x ~ a \ 
Then the number of inactive particles scales as Ni a ~ 
n(L) ~ £i-(i-a)/z acC ording to eq. (|17fl . which is in- 
deed a vanishing fraction of the active particles, which 
we have tacitly assumed. According to the latter expres- 
sion, if z + a < 1, Ni a — > as L — > 00, meaning that all 
particles are active. 

Setting a = in the scaling forms obtained for the 
current we recover the exact results known for the Sinai 
walk. 



V. ZERO AVERAGE BIAS 

If the control parameter is zero (5 = 0, the average tilt 
of the potential landscape, and consequently the sample- 
averaged current is zero even for finite systems. However, 
as a consequence of the fluctuations of randomness, the 
potential landscape of a single sample of size L is tilted, 
having a typical slope of L -1 / 2 , in accordance to central 
limit theorem. Therefore the quantity of interest is the 
magnitude of the current for a finite system. 

As the symmetric point is approached S — > 0, both 
the typical length and amplitude of Brownian excursions 
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are diverging as £ ~ 6~ 2 and U m ~ respectively. 
According to central limit theorem the typical amplitude 
is then of the order L 1 / 2 , consequently the magnitude of 
the current vanishes rapidly with the system size as 

-\n\J\~L x / 2 . (24) 

Since the largest excursion covers an 0(1) fraction of 
the landscape, its internal structure has to be investi- 
gated in order to get some insight into the properties 
of the density profile. We consider first a finite sample 
with $l = ^2i = ii^ n Pi ~ hi eft) = 0, where the potential 
is single-valued and the current is zero. In the station- 
ary state of this sample the particles practically occupy 
the N lowest-potential sites, and the clusters of occupied 
sites are represented by Brownian excursions. (These are 
not to be confused with barriers, which are excursions 
starting at extremal sites of the landscape.) The number 
of such excursions n in a system of size L is estimated as 
follows. For (5 = 0, the asymptotic distribution of their 
length is 

pi ~ r 3 / 2 , (25) 

thus the length of the largest one among n excursions, 
Imax, is typically of the order n 2 . On the other hand, n 
is related to L via n J^ max lp t dl ~ L, yielding n ~ L 1 ! 2 
and Imax ~ L. Thus, the potential landscape contains 
©(L 1 / 2 ) excursions and the largest one of them is macro- 
scopic. If the landscape is tilted, but L is large, one can 
still construct an approximate profile by cutting the land- 
scape by a straight line, the slope of which corresponds 
to the tilt of the landscape. Then sites below the line 
are filled with particles. Since the current as well as the 
typical tilt of samples vanish as L — > oo, the profile con- 
structed in this way is expected to be not much different 
from the true stationary profile. 

According to the above one expects that the distri- 
bution of interparticle distances has the asymptotics as 
given in (|25|) . with a finite cut-off 0(L). Indeed, nu- 
merical results are in satisfactory agreement with these 
predictions, see Fig. 

VI. RELATED MODELS 

In this section we shall apply the phenomenological 
description developed above to two generalizations of the 
disordered exclusion process. 

A. Particles of size d 

First, we consider a process, where the particles occupy 
d > 1 subsequent lattice sites whereas they take steps of 
unit length as before 3 ^. 

Obviously, the single-particle trapping time at a cer- 
tain barrier, r is independent of d. However, the time- 
scale for holes to overcome the same barrier which is now 
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FIG. 6: Distribution of interparticle distances calculated nu- 
merically with the bimodal distribution with c = 0.5 and 
r = 0.2. The number of samples is 5000. The slope of the 
straight line is —3/2. 

filled by particles of size d is proportional to r 1 ^, since 
the stepwidth of the holes is given by d, and the effec- 
tive potential difference they feel is divided by d. As- 
suming that there is a sharp front at the largest bar- 
rier, it must be located at the potential where the es- 
cape rate for particles r _1 and that for holes t _1 is re- 
lated as t _1 = dr^ 1 . Parameterizing the escape rate 
as t -1 = t 1 1 , where T\ is the single-particle trapping 
time, we have t 1 ^ l ^^ hd i _ ^ l / b _ p or \ ar g e Tl ^his 
yields b = 1 + d, consequently the current scales as 

J ~ L~ z ^ 1+d \ (26) 

where we made use of the scaling relation t± ~ L z . Fol- 
lowing the argumentions of the previous section, similar 
results can be derived here, as well. For instance, the cur- 
rent follows a Frechet distribution (see Fig0for numerical 
results) , and the number of inactive particles scales with 
the system size L as 

N. la - L d/{1+d) . (27) 

Finally, we mention, that the model with d — 2 is 
related to a PASEP in which the particles are linked, 
meaning that the number of empty sites between two 
neighboring particles can be at most one. It is easy to 
see that the motion of holes in this model follows the 
same rules as, say, the left half of d = 2 particles in the 
original model. In biological transport systems, molecu- 
lar motors are often attached to a rigid backbone, where 
the distance between motors is limited to a finite value 3 *. 
Exclusion processes with linked particles may be relevant 
for modelling cooperative behavior in these type of sys- 
tems. 



B. Exclusion process with multiple occupation 

Our second example is the disordered version of the 
generalized exclusion process 3 ^, where the number of par- 
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FIG. 7: Distribution of the current for the exclusion process 
with particles of size d = 2 (a) and d = 3 (b), computed for 
different system sizes. The number of particles is N = L/(2d). 
Bimodal randomness was used with c = 0.2, r = 0.3 where 
z = 0.869 (a) and c = 0.3, r = 0.3, where z = 1.420 (b), 
whereas the number of samples is 10 4 . The solid curves are 
the Frechet distributions given in 11611 . 



tides at a given site is limited by the site independent in- 
teger K > 1. The top most particles at a given site jumps 
to one of the neighboring sites with site-dependent rates 
Pi,qi, provided that the occupation at the target site is 
smaller than K. This model evidently interpolates be- 
tween the exclusion process (K — 1) and the zero-range 
process (K — N). Moreover, it can be roughly inter- 
preted as an exclusion process where the size of particles, 
d = 1/K, is smaller than one. 

In the stationary state of this generalized model we 
still expect a front at the largest barrier separating an 
almost fully occupied region (with K particles per site) 
from an almost empty region. Moreover, particle-hole 
symmetry holds for K > 1, as well, therefore half- filling 
must be realized at the largest barrier, and the current 
scales with the system size for any finite K as for K = 1 . 
Numerical results for the distribution of the current are 
presented in Fig. |H1 
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FIG. 8: Numerically calculated distributions of the current 
for the generalized exclusion process with K = 4, for different 
system sizes. The number of particles is iV = KL/2. Bimodal 
randomness was used with c = 0.2 and r — 0.5, where z = 0.5, 
whereas the number of samples is 10 4 . The solid curve is the 
Frechet distribution given in 1161 . 



One can also consider L-dependence in K, such as K = 
K L k , with < k < 1, leading to a scaling of the number 
of inactive particles scales according to Ni a ~ KL 1 / 2 ~ 
£i/2+*i resu j^ s f tfie previous section concerning 

the low-density limit can be generalized for this case in 
a straightforward way. 

We mention finally, that regarding the number of par- 
ticles at a given site as a distance variable, the present 
model can be mapped to an exclusion process with parti- 
clewise disorder, for any K , where the distance between 
neighboring particles can be at most K. 



VII. NON-STATIONARY PHENOMENA 

A. Coarsening 

In the following we analyse the approach to the sta- 
tionary state in an infinite system. When the system is 
started in a homogeneous configuration, it undergoes a 
coarsening process in the course of which the typical size 
of high-density and low-density regions, l(t), is growing 6 . 
The growth rate of these regions is determined by the 
current leading to the differential equation 



dl(t) 
dt 



\J{l{t))\- 



(28) 



where the current J (7(f)) is a function of the time- 
dependent length scale. By making use of the scaling 
relation of the current Ijl5(l which is valid for 5 ^ we 
get for large t 



Ut)~t 1/C , c = # + 1, (^o) 



(29) 



where £ is the dynamical exponent related to the coars- 
ening. We mention, that the same scaling relation has 
already been derived for the bimodal distribution, with 
z' instead of A, which is therefore valid only in the c — > 
limit. Using (|28|l and l|29l) . the asymptotic time depen- 
dence of the current is given by 



J(t) ~t~*Fv 



(6 jt 0). 



(30) 



When the symmetric point is approached, i.e. S — * 0, 
the exponent z given in l|12|) . and consequently the dy- 
namical exponent £ diverges as ( ~ z ~ l/(2<$)2&. Thus, 
strictly at the symmetric point (5 = 0, the length scale is 
growing slower than any power of t, a phenomenon called 
anomalous coarsening. Substituting the scaling relation 
of the current l|24|l into (|28|l . we get the asymptotic so- 
lution ^/ 2 e cil/2 



l(t) 



t, leading to 



t 



ln 'in^ 



(6 = 0) 



(31) 



which is indeed an anomalously (logarithmically) slow 
growth. Using H28JI we obtain for the asymptotic time- 
dependence of the magnitude of the current 



\m 



hxt 



(5 = 0) 



(32) 
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It is, however, circumstantial to measure the size of dense 
regions in numerical simulations since they are not con- 
nected. Instead, we considered the average displacement 
of particles (x(t)) as a function of time for S ^ 0, which is 
related to the current accor ding to = Com- 

paring this relation with (|28ll we see, that (x(t)} grows in 
time as the length scale l(t). Results of numerical simu- 
lations shown in Fig. El are in good agreement with the 
phenomenological predictions. 




4 8 12 16 

Int 



FIG. 9: Time dependence of the average displacement of a 
single particle measured for different values of 8 7^ 0. Bimodal 
randomness was used with c = | and r — -^=, ~, |, |, 

where z — 1,1,2,3,4, respectively. The size of the system 
is L = 8192 whereas TV = 4096 and the disorder average was 
performed over 200 samples. The straight lines have the slope 
^2 for each z. 

We mention, that the cut-off of the interparticle dis- 
tance distribution l(t) is growing slower than l(t). If 
2 < 1, the interparticle distances in the low-density re- 
gions of size l(t) are proportional to (l(t)) z / 2 , and thus 
grow as J(t) ~ t z /( 2+z \ If 1 < z < oo, the interparticle 
distances in the low-density regions are proportional to 
{lit)) 1 / 2 leading to l{t) ~ tVP+f). 



B. Invasion 

We consider in this section an open semi-infinite lattice 
with entrance rate a = 1, furthermore, it is assumed 
that the lattice is initially empty. Our interest is in the 
invading dynamics, i.e. the dynamics of the first particle 
that enters the system and the motion of the bulk of 
particles. This question is of practical relevance, e.g. one 
considers a fluid penetrating a porous medium. 

First, we study how the position of the leading particle 
evolves in time. Apparently, it must not travel slower 
than a single walker since the jumps to the right are 
never hindered, contrary to jumps to the left. Thus, for 
z < 1 the first particle advances with a constant velocity. 
For z > 1, we consider a finite system of size L, and 
determine the finite-size scaling of the characteristic time 



t p needed for the first particle to reach site L, as follows. 
The leading contribution to r p comes from the trapping 
time at the largest barrier and at subsequent barriers. 
Therefore we assume that the first particle has already 
arrived at the largest barrier of the system. Without the 
other particles, the typical time to overcome the barrier 
would be 7~i ~ L z . However, during this period other 
particles are arriving at the barrier with rate 0{L~ z l 2 ) 
determined by the current at the half-filled largest barrier 
between the entrance site and leading particle. The size 
of the largest barrier is at most ~ InL, so it takes at 
most a period of t ~ L z / 2 lnL until half-filling of the 
largest barrier is achieved. Once half-filling is realized the 
trapping time for the leading particle is only 0(L Z / 2 ), so 
the total trapping time at the largest barrier is at most 
0{L Z / 2 InL). 

After leaving the largest barrier, the first particle ad- 
vances freely until it arrives at a barrier with trapping 
time larger then a/tJ, where it is assisted in passing the 
barrier again by other particles. The typical time-scale to 
fill the barrier is determined by the current and therefore 
given by t ~ L z l 2 . Using the assistance of the following 
particles the trapping time of the first particle at the bar- 
rier is given by 0(L z l 2 ). Since there are 0{L X / 2 ) barriers 
with trapping time r > ^Jt\ in a sample of size L, the 
characteristic time for getting through the system scales 
in leading order as r p ~ i 2 / 2 + 1 / 2 . Therefore in a semi- 
infinite system, the position of the leading particle £ p is 
expected to scale as 

£ p ~t 1 /'*, with Zp = ^p-, (33) 

for z > 1, whereas z p — 1 for z < 1. 

In order to test the validity of our approach, we have 
checked the above results numerically by measuring the 
distribution of times the first particle needs to traverse 
a finite system. The results are in good agreement with 
the scaling relation (|33|l . as shown in Fig. For zero 
average bias, the exponent z p is formally infinite and the 
propagation of the first particle is anomalously slow. 
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FIG. 10: Distribution of invasion times calculated for differ- 
ent system sizes. Bimodal randomness was used with c = 0.3 
and r = 0.2, where z p ~ 1.45. The number of samples is 10 4 . 

Compared to the first particle, the bulk of particles 
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moves slower, as we shall show below. For large times, 
t> 1, a dense particle cluster is present in the system, 
which extends from the entrance site to a barrier which is 
located at and which must be the largest barrier of 
the interval [0, £/(£)]. Since this barrier is half-filled, the 
current it maintains is typically in the order of £,f(t)~ z ^ 2 . 
Through this (source) barrier particles are constantly 
moving to the next larger barrier to its right. The es- 
cape rate at this 2nd barrier is typically (2£/(i))~ z,/2 , 
which is a finite fraction of the filling current, therefore 
the domain behind it is being filled up by particles with 
a rate ~ J ~ £,}{t)~ z ^ 2 - (Those particles which escape 
at the 2nd barrier meanwhile, rush forward, to the third 
and subsequent larger barriers spreading in the domain 
behind the leading particle.) After certain time the do- 
main between the source and the 2nd barrier is com- 
pletely filled, and the new position of the front is now at 
the 2nd barrier. Thus the growth rate of the dense clus- 
ter is proportional to the current, ~ J ~ £f{t)~ z ^ 2 
yielding that the typical position of the front scales with 
the time as 

* 1/C , (34) 

with a larger dynamical exponent than that of the lead- 
ing particle. In order to check this numerically we have 
measured the number of particles N(t) as a function of 
time, which is expected to grow as £/. Results of numer- 
ical simulations shown in Fig 1111 are in good agreement 
with QOJ. 
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FIG. 11: Time dependence of the average number of particles 
of measured for different values of 5 7^ 0. Bimodal random- 
ness was used with c = ^ and r = |, j, g where z = 1, 2, 3, 
respectively. The size of the system is L = 8192 and the dis- 
order average was performed over 200 samples. The straight 
lines have the slope for each z. 



VIII. DISCUSSION 

In this paper we have studied the partially asymmet- 
ric exclusion process with lattice disorder and we have 



obtained several conjecturedly exact results both for the 
stationary state and for related non-stationary problems. 
In particular for the biased system we have shown a 
phase separation phenomenon where a domain wall tak- 
ing place in the middle of the largest barrier separates 
a low- and a high-density phase. The stationary current 
in a large finite system goes to zero as a power-law, and 
the corresponding exponent is just the half of the dy- 
namical exponent of a single random walker. We have 
also classified the particles, which typically belong to a 
macroscopic condensate. This condensate, however, is 
cut up by O^L 1 / 2 ) vacant regions of finite size, which 
are located at the subleading barriers of the potential 
landscape. These vacant regions, through the particle- 
hole symmetry correspond to filled deep valleys in the 
low-density phase, which contain the so called inactive 
particles. Finally, the third type of particles are the ac- 
tive ones, which actually carry the current in the sys- 
tem. We have also studied the properties of the system 
in non-stationary processes, such as in coarsening or in 
the invasion, and have shown that the scaling exponents 
are related to the dynamical exponent of a single random 
walker. 

Comparing the effect of different types of disorder on 
the stationary behavior of the PASEP we can notice some 
analogies but also several differences. Both for particle 
and lattice disorder the stationary current is vanishing 
as a power-law of the size of the lattice, and the crit- 
ical exponents are simple related to each other. Also 
the phase-separation phenomena is common for the two 
problems. The system with particle disorder can be 
treated by renormalizing the particles^ arriving to an ef- 
fective model containing only a single particle. Renor- 
malization for lattice disorder has only been performed 
for a single particle 2 ^ leading to a motion of the particle 
in a renormalized landscape. For many particles it is ex- 
pected that both the landscape and the particles should 
be renormalized and the resulting model is still a many 
particle problem. The actual construction of the renor- 
malization procedure is the purpose of future research. 
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